import mne
import matplotlib.pyplot as plt
%matplotlib inline
from metadata import *
print(subject_list)
print(EOG_channels)
['001', '002', '003', '004', '005', '006', '007', '008', '009', '010', '011', '012', '015', '016', '021', '022', '024', '025', '027', '029', '030', '031', '032', '033'] ['EXG3', 'EXG4', 'EXG5'] C:\Users\DELL\AppData\Local\Programs\Python\Python37\lib\site-packages\mne\utils\docs.py:830: DeprecationWarning: Function read_montage is deprecated; ``read_montage`` is deprecated and will be removed in v0.20. Please use ``read_dig_fif``, ``read_dig_egi``, ``read_custom_montage``, or ``read_dig_captrack`` to read a digitization based on your needs instead; or ``make_standard_montage`` to create ``DigMontage`` based on template; or ``make_dig_montage`` to create a ``DigMontage`` out of np.arrays warnings.warn(msg, category=DeprecationWarning) C:\Users\DELL\AppData\Local\Programs\Python\Python37\lib\site-packages\mne\utils\docs.py:813: DeprecationWarning: Class Montage is deprecated; Montage class is deprecated and will be removed in v0.20. Please use DigMontage instead. warnings.warn(msg, category=DeprecationWarning)
# HELPTER FUNCTIONS
def read_data(subject):
print('# Reading Raw and ICA data...')
raw = mne.io.read_raw_fif(filtered_files[subject], preload=True)
ica = mne.preprocessing.read_ica(ica_files[subject])
ica.exclude = [] # default to empty list
return raw, ica
def show_ics(subject):
%matplotlib qt
raw, ica = read_data(subject)
ica.plot_sources(raw)
# read more on https://mne.tools/dev/generated/mne.preprocessing.ICA.html#mne.preprocessing.ICA.find_bads_eog
def auto_detect_eog(subject, raw, ica, EOG):
print('# Finding ICs that match patterns in EOG channels...')
exclude = set()
for channel in EOG:
if (channel in raw.info['bads']): # skip bad channels
continue
eog_indices, eog_scores = ica.find_bads_eog(
raw,
ch_name = channel
)
ica.plot_scores(eog_scores)
exclude.update(eog_indices)
print('# Exclude EOG:', exclude)
return exclude
# read more on https://mne.tools/dev/generated/mne.preprocessing.ICA.html#mne.preprocessing.ICA.find_bads_ecg
def auto_detect_ecg(subject, raw, ica, ECG):
print('# Finding ICs that match patterns in ECG channels...')
exclude = set()
for channel in ECG:
if (channel in raw.info['bads']): # skip bad channels
continue
ecg_indices, ecg_scores = ica.find_bads_ecg(
raw,
ch_name = channel
)
ica.plot_scores(ecg_scores)
exclude.update(ecg_indices)
print('# Exclude ECG:', exclude)
return exclude
def reject(subject, EOG, ECG, noise):
%matplotlib inline
print('##### Processing subject', subject, '...')
raw, ica = read_data(subject)
print('# Running auto detection...')
exclude_eog = auto_detect_eog(subject, raw, ica, EOG)
exclude_ecg = auto_detect_ecg(subject, raw, ica, ECG)
exclude_auto = exclude_eog.union(exclude_eog)
print('# Auto Exclude:', exclude_auto)
exclude_list = list(exclude_auto.union(set(noise))) # automatically detected + manually included
print('# Exclude list:', exclude_list)
ica.exclude = exclude_list
print('# Visualizing noise reduction...')
picks = EEG_channels
for bad in raw.info['bads']:
if (bad in picks):
picks.remove(bad)
ica.plot_overlay(raw, picks = picks)
print("# Saving the updated data...")
raw_after_ica = ica.apply(raw, exclude = ica.exclude)
raw_after_ica.save(ica_after_files[subject], overwrite = True)
print('Done!')
'''
FOR EACH SUBJECT, DROP BAD COMPONENTS.
subejct: subject_id
EOG_list: already specified by EOG_channels in "metadata" file: ['EXG3', 'EXG4', 'EXG5']
ECG_list: default to be empty
noise_list: manually specified by inspecting visualizations of independent components
finally just run: reject(subject, EOG_channels, ECG_list, noise_list)
'''
'\nFOR EACH SUBJECT, DROP BAD COMPONENTS. \n\nsubejct: subject_id\nEOG_list: already specified by EOG_channels in "metadata" file: [\'EXG3\', \'EXG4\', \'EXG5\']\nECG_list: default to be empty\nnoise_list: manually specified by inspecting visualizations of independent components\n\nfinally just run: reject(subject, EOG_channels, ECG_list, noise_list)\n'
subject = '001'
ECG_list = []
noise_list = [1, 5, 7, 8, 9, 16, 35, 36] # eye: 1, 5, 7; muscle: 8, 9; noise: 35, 36
reject(subject, EOG_channels, ECG_list, noise_list)
##### Processing subject 001 ...
# Reading Raw and ICA data...
Opening raw data file F:/HoltLab/data/filtered_raw/001_filtered_raw.fif...
Isotrak not found
Range : 0 ... 638591 = 0.000 ... 4988.992 secs
Ready.
Reading 0 ... 638591 = 0.000 ... 4988.992 secs...
Reading F:/HoltLab/data/ica/001_ica.fif ...
Now restoring ICA solution ...
Ready.
# Running auto detection...
# Finding ICs that match patterns in EOG channels...
Using channel EXG3 as EOG channel
... filtering ICA sources
Setting up band-pass filter from 1 - 10 Hz
FIR filter parameters
---------------------
Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter:
- Windowed frequency-domain design (firwin2) method
- Hann window
- Lower passband edge: 1.00
- Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz)
- Upper passband edge: 10.00 Hz
- Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz)
- Filter length: 2048 samples (16.000 sec)
... filtering target
Setting up band-pass filter from 1 - 10 Hz
FIR filter parameters
---------------------
Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter:
- Windowed frequency-domain design (firwin2) method
- Hann window
- Lower passband edge: 1.00
- Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz)
- Upper passband edge: 10.00 Hz
- Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz)
- Filter length: 2048 samples (16.000 sec)
Using channel EXG4 as EOG channel ... filtering ICA sources Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec) ... filtering target Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec)
Using channel EXG5 as EOG channel ... filtering ICA sources Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec) ... filtering target Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec)
# Exclude EOG: {0, 4, 36}
# Finding ICs that match patterns in ECG channels...
# Exclude ECG: set()
# Auto Exclude: {0, 4, 36}
# Exclude list: [0, 1, 35, 4, 36, 5, 7, 8, 9, 16]
# Visualizing noise reduction...
Transforming to ICA space (37 components)
Zeroing out 10 ICA components
# Saving the updated data... Transforming to ICA space (37 components) Zeroing out 10 ICA components Overwriting existing file. Writing F:\HoltLab\data\raw_after_ica\001_after_ica.fif Closing F:\HoltLab\data\raw_after_ica\001_after_ica.fif [done] Done!
subject = '002'
ECG_list = []
noise_list = [1, 4, 8, 9, 28, 36] # eye: 1; muscle: 4, 8, 9, 28; noise: 36
reject(subject, EOG_channels, ECG_list, noise_list)
##### Processing subject 002 ...
# Reading Raw and ICA data...
Opening raw data file F:/HoltLab/data/filtered_raw/002_filtered_raw.fif...
Isotrak not found
Range : 0 ... 633343 = 0.000 ... 4947.992 secs
Ready.
Reading 0 ... 633343 = 0.000 ... 4947.992 secs...
Reading F:/HoltLab/data/ica/002_ica.fif ...
Now restoring ICA solution ...
Ready.
# Running auto detection...
# Finding ICs that match patterns in EOG channels...
Using channel EXG3 as EOG channel
... filtering ICA sources
Setting up band-pass filter from 1 - 10 Hz
FIR filter parameters
---------------------
Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter:
- Windowed frequency-domain design (firwin2) method
- Hann window
- Lower passband edge: 1.00
- Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz)
- Upper passband edge: 10.00 Hz
- Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz)
- Filter length: 2048 samples (16.000 sec)
... filtering target
Setting up band-pass filter from 1 - 10 Hz
FIR filter parameters
---------------------
Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter:
- Windowed frequency-domain design (firwin2) method
- Hann window
- Lower passband edge: 1.00
- Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz)
- Upper passband edge: 10.00 Hz
- Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz)
- Filter length: 2048 samples (16.000 sec)
Using channel EXG4 as EOG channel ... filtering ICA sources Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec) ... filtering target Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec)
Using channel EXG5 as EOG channel ... filtering ICA sources Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec) ... filtering target Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec)
# Exclude EOG: {0, 36, 6}
# Finding ICs that match patterns in ECG channels...
# Exclude ECG: set()
# Auto Exclude: {0, 36, 6}
# Exclude list: [0, 1, 36, 4, 6, 8, 9, 28]
# Visualizing noise reduction...
Transforming to ICA space (37 components)
Zeroing out 8 ICA components
# Saving the updated data... Transforming to ICA space (37 components) Zeroing out 8 ICA components Overwriting existing file. Writing F:\HoltLab\data\raw_after_ica\002_after_ica.fif Closing F:\HoltLab\data\raw_after_ica\002_after_ica.fif [done] Done!
subject = '003'
ECG_list = []
noise_list = [6, 7, 12, 13, 17, 22, 23, 26, 28, 29, 36] # eye: 6, 17; muscle: 7, 12, 13, 22, 36; noise: 23, 26, 28, 29
reject(subject, EOG_channels, ECG_list, noise_list)
##### Processing subject 003 ...
# Reading Raw and ICA data...
Opening raw data file F:/HoltLab/data/filtered_raw/003_filtered_raw.fif...
Isotrak not found
Range : 0 ... 615295 = 0.000 ... 4806.992 secs
Ready.
Reading 0 ... 615295 = 0.000 ... 4806.992 secs...
Reading F:/HoltLab/data/ica/003_ica.fif ...
Now restoring ICA solution ...
Ready.
# Running auto detection...
# Finding ICs that match patterns in EOG channels...
Using channel EXG3 as EOG channel
... filtering ICA sources
Setting up band-pass filter from 1 - 10 Hz
FIR filter parameters
---------------------
Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter:
- Windowed frequency-domain design (firwin2) method
- Hann window
- Lower passband edge: 1.00
- Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz)
- Upper passband edge: 10.00 Hz
- Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz)
- Filter length: 2048 samples (16.000 sec)
... filtering target
Setting up band-pass filter from 1 - 10 Hz
FIR filter parameters
---------------------
Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter:
- Windowed frequency-domain design (firwin2) method
- Hann window
- Lower passband edge: 1.00
- Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz)
- Upper passband edge: 10.00 Hz
- Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz)
- Filter length: 2048 samples (16.000 sec)
Using channel EXG4 as EOG channel ... filtering ICA sources Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec) ... filtering target Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec)
Using channel EXG5 as EOG channel ... filtering ICA sources Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec) ... filtering target Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec)
# Exclude EOG: {0, 4}
# Finding ICs that match patterns in ECG channels...
# Exclude ECG: set()
# Auto Exclude: {0, 4}
# Exclude list: [0, 4, 36, 6, 7, 12, 13, 17, 22, 23, 26, 28, 29]
# Visualizing noise reduction...
Transforming to ICA space (37 components)
Zeroing out 13 ICA components
# Saving the updated data... Transforming to ICA space (37 components) Zeroing out 13 ICA components Overwriting existing file. Writing F:\HoltLab\data\raw_after_ica\003_after_ica.fif Closing F:\HoltLab\data\raw_after_ica\003_after_ica.fif [done] Done!
subject = '004'
ECG_list = []
noise_list = [7, 8, 12, 20, 21, 23, 35] # eye: 7, 12, 20, 35; noise: 8, 21, 23
reject(subject, EOG_channels, ECG_list, noise_list)
##### Processing subject 004 ...
# Reading Raw and ICA data...
Opening raw data file F:/HoltLab/data/filtered_raw/004_filtered_raw.fif...
Isotrak not found
Range : 0 ... 643327 = 0.000 ... 5025.992 secs
Ready.
Reading 0 ... 643327 = 0.000 ... 5025.992 secs...
Reading F:/HoltLab/data/ica/004_ica.fif ...
Now restoring ICA solution ...
Ready.
# Running auto detection...
# Finding ICs that match patterns in EOG channels...
Using channel EXG3 as EOG channel
... filtering ICA sources
Setting up band-pass filter from 1 - 10 Hz
FIR filter parameters
---------------------
Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter:
- Windowed frequency-domain design (firwin2) method
- Hann window
- Lower passband edge: 1.00
- Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz)
- Upper passband edge: 10.00 Hz
- Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz)
- Filter length: 2048 samples (16.000 sec)
... filtering target
Setting up band-pass filter from 1 - 10 Hz
FIR filter parameters
---------------------
Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter:
- Windowed frequency-domain design (firwin2) method
- Hann window
- Lower passband edge: 1.00
- Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz)
- Upper passband edge: 10.00 Hz
- Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz)
- Filter length: 2048 samples (16.000 sec)
Using channel EXG4 as EOG channel ... filtering ICA sources Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec) ... filtering target Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec)
Using channel EXG5 as EOG channel ... filtering ICA sources Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec) ... filtering target Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec)
# Exclude EOG: {0, 1}
# Finding ICs that match patterns in ECG channels...
# Exclude ECG: set()
# Auto Exclude: {0, 1}
# Exclude list: [0, 1, 35, 7, 8, 12, 20, 21, 23]
# Visualizing noise reduction...
Transforming to ICA space (36 components)
Zeroing out 9 ICA components
# Saving the updated data... Transforming to ICA space (36 components) Zeroing out 9 ICA components Overwriting existing file. Writing F:\HoltLab\data\raw_after_ica\004_after_ica.fif Closing F:\HoltLab\data\raw_after_ica\004_after_ica.fif [done] Done!
subject = '005'
ECG_list = []
noise_list = [0, 1, 3, 9, 22, 30, 31, 34, 36] # muscle: 22, 31, 34, 36; noise: 0, 1, 3, 9, 30
reject(subject, EOG_channels, ECG_list, noise_list)
##### Processing subject 005 ...
# Reading Raw and ICA data...
Opening raw data file F:/HoltLab/data/filtered_raw/005_filtered_raw.fif...
Isotrak not found
Range : 0 ... 627711 = 0.000 ... 4903.992 secs
Ready.
Reading 0 ... 627711 = 0.000 ... 4903.992 secs...
Reading F:/HoltLab/data/ica/005_ica.fif ...
Now restoring ICA solution ...
Ready.
# Running auto detection...
# Finding ICs that match patterns in EOG channels...
Using channel EXG3 as EOG channel
... filtering ICA sources
Setting up band-pass filter from 1 - 10 Hz
FIR filter parameters
---------------------
Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter:
- Windowed frequency-domain design (firwin2) method
- Hann window
- Lower passband edge: 1.00
- Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz)
- Upper passband edge: 10.00 Hz
- Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz)
- Filter length: 2048 samples (16.000 sec)
... filtering target
Setting up band-pass filter from 1 - 10 Hz
FIR filter parameters
---------------------
Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter:
- Windowed frequency-domain design (firwin2) method
- Hann window
- Lower passband edge: 1.00
- Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz)
- Upper passband edge: 10.00 Hz
- Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz)
- Filter length: 2048 samples (16.000 sec)
Using channel EXG4 as EOG channel ... filtering ICA sources Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec) ... filtering target Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec)
Using channel EXG5 as EOG channel ... filtering ICA sources Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec) ... filtering target Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec)
# Exclude EOG: {2, 4, 6}
# Finding ICs that match patterns in ECG channels...
# Exclude ECG: set()
# Auto Exclude: {2, 4, 6}
# Exclude list: [0, 1, 2, 34, 4, 3, 6, 36, 9, 22, 30, 31]
# Visualizing noise reduction...
Transforming to ICA space (37 components)
Zeroing out 12 ICA components
# Saving the updated data... Transforming to ICA space (37 components) Zeroing out 12 ICA components Overwriting existing file. Writing F:\HoltLab\data\raw_after_ica\005_after_ica.fif Closing F:\HoltLab\data\raw_after_ica\005_after_ica.fif [done] Done!
subject = '006'
ECG_list = []
noise_list = [2, 3, 4, 5, 10, 18, 19, 20, 24, 25, 34] # eye: 2, 4, 20, 24; muscle: 10, 18, 19, 25, 34; noise: 3, 5
reject(subject, EOG_channels, ECG_list, noise_list)
##### Processing subject 006 ...
# Reading Raw and ICA data...
Opening raw data file F:/HoltLab/data/filtered_raw/006_filtered_raw.fif...
Isotrak not found
Range : 0 ... 843647 = 0.000 ... 6590.992 secs
Ready.
Reading 0 ... 843647 = 0.000 ... 6590.992 secs...
Reading F:/HoltLab/data/ica/006_ica.fif ...
Now restoring ICA solution ...
Ready.
# Running auto detection...
# Finding ICs that match patterns in EOG channels...
Using channel EXG3 as EOG channel
... filtering ICA sources
Setting up band-pass filter from 1 - 10 Hz
FIR filter parameters
---------------------
Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter:
- Windowed frequency-domain design (firwin2) method
- Hann window
- Lower passband edge: 1.00
- Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz)
- Upper passband edge: 10.00 Hz
- Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz)
- Filter length: 2048 samples (16.000 sec)
... filtering target
Setting up band-pass filter from 1 - 10 Hz
FIR filter parameters
---------------------
Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter:
- Windowed frequency-domain design (firwin2) method
- Hann window
- Lower passband edge: 1.00
- Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz)
- Upper passband edge: 10.00 Hz
- Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz)
- Filter length: 2048 samples (16.000 sec)
Using channel EXG4 as EOG channel ... filtering ICA sources Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec) ... filtering target Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec)
Using channel EXG5 as EOG channel ... filtering ICA sources Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec) ... filtering target Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec)
# Exclude EOG: {0, 1}
# Finding ICs that match patterns in ECG channels...
# Exclude ECG: set()
# Auto Exclude: {0, 1}
# Exclude list: [0, 1, 2, 3, 4, 5, 34, 10, 18, 19, 20, 24, 25]
# Visualizing noise reduction...
Transforming to ICA space (35 components)
Zeroing out 13 ICA components
# Saving the updated data... Transforming to ICA space (35 components) Zeroing out 13 ICA components Overwriting existing file. Writing F:\HoltLab\data\raw_after_ica\006_after_ica.fif Closing F:\HoltLab\data\raw_after_ica\006_after_ica.fif [done] Done!
subject = '007'
ECG_list = []
noise_list = [12, 16, 19, 22, 28] # eye: 16, 19, 28; muscle: 12; noise: 22
reject(subject, EOG_channels, ECG_list, noise_list)
##### Processing subject 007 ...
# Reading Raw and ICA data...
Opening raw data file F:/HoltLab/data/filtered_raw/007_filtered_raw.fif...
Isotrak not found
Range : 0 ... 779263 = 0.000 ... 6087.992 secs
Ready.
Reading 0 ... 779263 = 0.000 ... 6087.992 secs...
Reading F:/HoltLab/data/ica/007_ica.fif ...
Now restoring ICA solution ...
Ready.
# Running auto detection...
# Finding ICs that match patterns in EOG channels...
Using channel EXG3 as EOG channel
... filtering ICA sources
Setting up band-pass filter from 1 - 10 Hz
FIR filter parameters
---------------------
Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter:
- Windowed frequency-domain design (firwin2) method
- Hann window
- Lower passband edge: 1.00
- Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz)
- Upper passband edge: 10.00 Hz
- Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz)
- Filter length: 2048 samples (16.000 sec)
... filtering target
Setting up band-pass filter from 1 - 10 Hz
FIR filter parameters
---------------------
Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter:
- Windowed frequency-domain design (firwin2) method
- Hann window
- Lower passband edge: 1.00
- Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz)
- Upper passband edge: 10.00 Hz
- Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz)
- Filter length: 2048 samples (16.000 sec)
Using channel EXG4 as EOG channel ... filtering ICA sources Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec) ... filtering target Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec)
Using channel EXG5 as EOG channel ... filtering ICA sources Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec) ... filtering target Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec)
# Exclude EOG: {0, 1, 2, 4}
# Finding ICs that match patterns in ECG channels...
# Exclude ECG: set()
# Auto Exclude: {0, 1, 2, 4}
# Exclude list: [0, 1, 2, 4, 12, 16, 19, 22, 28]
# Visualizing noise reduction...
Transforming to ICA space (37 components)
Zeroing out 9 ICA components
# Saving the updated data... Transforming to ICA space (37 components) Zeroing out 9 ICA components Overwriting existing file. Writing F:\HoltLab\data\raw_after_ica\007_after_ica.fif Closing F:\HoltLab\data\raw_after_ica\007_after_ica.fif [done] Done!
subject = '008'
ECG_list = []
noise_list = [6, 19, 20, 26, 31, 32] # eye: 6, 26; muscle: 20, 31, 32; noise: 19
reject(subject, EOG_channels, ECG_list, noise_list)
##### Processing subject 008 ...
# Reading Raw and ICA data...
Opening raw data file F:/HoltLab/data/filtered_raw/008_filtered_raw.fif...
Isotrak not found
Range : 0 ... 629759 = 0.000 ... 4919.992 secs
Ready.
Reading 0 ... 629759 = 0.000 ... 4919.992 secs...
Reading F:/HoltLab/data/ica/008_ica.fif ...
Now restoring ICA solution ...
Ready.
# Running auto detection...
# Finding ICs that match patterns in EOG channels...
Using channel EXG3 as EOG channel
... filtering ICA sources
Setting up band-pass filter from 1 - 10 Hz
FIR filter parameters
---------------------
Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter:
- Windowed frequency-domain design (firwin2) method
- Hann window
- Lower passband edge: 1.00
- Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz)
- Upper passband edge: 10.00 Hz
- Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz)
- Filter length: 2048 samples (16.000 sec)
... filtering target
Setting up band-pass filter from 1 - 10 Hz
FIR filter parameters
---------------------
Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter:
- Windowed frequency-domain design (firwin2) method
- Hann window
- Lower passband edge: 1.00
- Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz)
- Upper passband edge: 10.00 Hz
- Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz)
- Filter length: 2048 samples (16.000 sec)
Using channel EXG4 as EOG channel ... filtering ICA sources Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec) ... filtering target Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec)
Using channel EXG5 as EOG channel ... filtering ICA sources Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec) ... filtering target Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec)
# Exclude EOG: {2, 10, 27}
# Finding ICs that match patterns in ECG channels...
# Exclude ECG: set()
# Auto Exclude: {2, 10, 27}
# Exclude list: [32, 2, 6, 10, 19, 20, 26, 27, 31]
# Visualizing noise reduction...
Transforming to ICA space (37 components)
Zeroing out 9 ICA components
# Saving the updated data... Transforming to ICA space (37 components) Zeroing out 9 ICA components Overwriting existing file. Writing F:\HoltLab\data\raw_after_ica\008_after_ica.fif Closing F:\HoltLab\data\raw_after_ica\008_after_ica.fif [done] Done!
subject = '009'
ECG_list = []
noise_list = [0, 2, 3, 4, 8, 9, 11, 19, 24] # eye: 0, 4; muscle: 2, 3, 8, 9, 11, 24; noise: 19
reject(subject, EOG_channels, ECG_list, noise_list)
##### Processing subject 009 ...
# Reading Raw and ICA data...
Opening raw data file F:/HoltLab/data/filtered_raw/009_filtered_raw.fif...
Isotrak not found
Range : 0 ... 710911 = 0.000 ... 5553.992 secs
Ready.
Reading 0 ... 710911 = 0.000 ... 5553.992 secs...
Reading F:/HoltLab/data/ica/009_ica.fif ...
Now restoring ICA solution ...
Ready.
# Running auto detection...
# Finding ICs that match patterns in EOG channels...
Using channel EXG3 as EOG channel
... filtering ICA sources
Setting up band-pass filter from 1 - 10 Hz
FIR filter parameters
---------------------
Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter:
- Windowed frequency-domain design (firwin2) method
- Hann window
- Lower passband edge: 1.00
- Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz)
- Upper passband edge: 10.00 Hz
- Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz)
- Filter length: 2048 samples (16.000 sec)
... filtering target
Setting up band-pass filter from 1 - 10 Hz
FIR filter parameters
---------------------
Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter:
- Windowed frequency-domain design (firwin2) method
- Hann window
- Lower passband edge: 1.00
- Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz)
- Upper passband edge: 10.00 Hz
- Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz)
- Filter length: 2048 samples (16.000 sec)
Using channel EXG4 as EOG channel ... filtering ICA sources Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec) ... filtering target Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec)
Using channel EXG5 as EOG channel ... filtering ICA sources Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec) ... filtering target Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec)
# Exclude EOG: {1, 5}
# Finding ICs that match patterns in ECG channels...
# Exclude ECG: set()
# Auto Exclude: {1, 5}
# Exclude list: [0, 1, 2, 3, 4, 5, 8, 9, 11, 19, 24]
# Visualizing noise reduction...
Transforming to ICA space (37 components)
Zeroing out 11 ICA components
# Saving the updated data... Transforming to ICA space (37 components) Zeroing out 11 ICA components Overwriting existing file. Writing F:\HoltLab\data\raw_after_ica\009_after_ica.fif Closing F:\HoltLab\data\raw_after_ica\009_after_ica.fif [done] Done!
subject = '010'
ECG_list = []
noise_list = [0, 3, 9, 11, 12, 13, 14, 17, 20, 21, 22, 23] # eye: 9, 11, 13, 20, 21; muscle: 12, 14, 17, 23; noise: 0, 3, 22
reject(subject, EOG_channels, ECG_list, noise_list)
##### Processing subject 010 ...
# Reading Raw and ICA data...
Opening raw data file F:/HoltLab/data/filtered_raw/010_filtered_raw.fif...
Isotrak not found
Range : 0 ... 679039 = 0.000 ... 5304.992 secs
Ready.
Reading 0 ... 679039 = 0.000 ... 5304.992 secs...
Reading F:/HoltLab/data/ica/010_ica.fif ...
Now restoring ICA solution ...
Ready.
# Running auto detection...
# Finding ICs that match patterns in EOG channels...
Using channel EXG3 as EOG channel
... filtering ICA sources
Setting up band-pass filter from 1 - 10 Hz
FIR filter parameters
---------------------
Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter:
- Windowed frequency-domain design (firwin2) method
- Hann window
- Lower passband edge: 1.00
- Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz)
- Upper passband edge: 10.00 Hz
- Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz)
- Filter length: 2048 samples (16.000 sec)
... filtering target
Setting up band-pass filter from 1 - 10 Hz
FIR filter parameters
---------------------
Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter:
- Windowed frequency-domain design (firwin2) method
- Hann window
- Lower passband edge: 1.00
- Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz)
- Upper passband edge: 10.00 Hz
- Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz)
- Filter length: 2048 samples (16.000 sec)
Using channel EXG4 as EOG channel ... filtering ICA sources Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec) ... filtering target Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec)
Using channel EXG5 as EOG channel ... filtering ICA sources Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec) ... filtering target Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec)
# Exclude EOG: {1, 2, 4, 5, 6}
# Finding ICs that match patterns in ECG channels...
# Exclude ECG: set()
# Auto Exclude: {1, 2, 4, 5, 6}
# Exclude list: [0, 1, 2, 3, 4, 5, 6, 9, 11, 12, 13, 14, 17, 20, 21, 22, 23]
# Visualizing noise reduction...
Transforming to ICA space (36 components)
Zeroing out 17 ICA components
# Saving the updated data... Transforming to ICA space (36 components) Zeroing out 17 ICA components Overwriting existing file. Writing F:\HoltLab\data\raw_after_ica\010_after_ica.fif Closing F:\HoltLab\data\raw_after_ica\010_after_ica.fif [done] Done!
subject = '011'
ECG_list = []
noise_list = [1, 11, 16, 19, 20, 24, 27, 33] # eye: 1; muscle: 11, 16, 19, 20, 24, 27, 33
reject(subject, EOG_channels, ECG_list, noise_list)
##### Processing subject 011 ...
# Reading Raw and ICA data...
Opening raw data file F:/HoltLab/data/filtered_raw/011_filtered_raw.fif...
Isotrak not found
Range : 0 ... 686975 = 0.000 ... 5366.992 secs
Ready.
Reading 0 ... 686975 = 0.000 ... 5366.992 secs...
Reading F:/HoltLab/data/ica/011_ica.fif ...
Now restoring ICA solution ...
Ready.
# Running auto detection...
# Finding ICs that match patterns in EOG channels...
Using channel EXG3 as EOG channel
... filtering ICA sources
Setting up band-pass filter from 1 - 10 Hz
FIR filter parameters
---------------------
Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter:
- Windowed frequency-domain design (firwin2) method
- Hann window
- Lower passband edge: 1.00
- Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz)
- Upper passband edge: 10.00 Hz
- Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz)
- Filter length: 2048 samples (16.000 sec)
... filtering target
Setting up band-pass filter from 1 - 10 Hz
FIR filter parameters
---------------------
Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter:
- Windowed frequency-domain design (firwin2) method
- Hann window
- Lower passband edge: 1.00
- Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz)
- Upper passband edge: 10.00 Hz
- Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz)
- Filter length: 2048 samples (16.000 sec)
Using channel EXG4 as EOG channel ... filtering ICA sources Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec) ... filtering target Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec)
Using channel EXG5 as EOG channel ... filtering ICA sources Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec) ... filtering target Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec)
# Exclude EOG: {9, 3, 6}
# Finding ICs that match patterns in ECG channels...
# Exclude ECG: set()
# Auto Exclude: {9, 3, 6}
# Exclude list: [1, 33, 3, 6, 9, 11, 16, 19, 20, 24, 27]
# Visualizing noise reduction...
Transforming to ICA space (34 components)
Zeroing out 11 ICA components
# Saving the updated data... Transforming to ICA space (34 components) Zeroing out 11 ICA components Overwriting existing file. Writing F:\HoltLab\data\raw_after_ica\011_after_ica.fif Closing F:\HoltLab\data\raw_after_ica\011_after_ica.fif [done] Done!
subject = '012'
ECG_list = []
noise_list = [3, 4, 7, 10, 25] # eye: 3, 4, 7, 25; muscle: 10
reject(subject, EOG_channels, ECG_list, noise_list)
##### Processing subject 012 ...
# Reading Raw and ICA data...
Opening raw data file F:/HoltLab/data/filtered_raw/012_filtered_raw.fif...
Isotrak not found
Range : 0 ... 628863 = 0.000 ... 4912.992 secs
Ready.
Reading 0 ... 628863 = 0.000 ... 4912.992 secs...
Reading F:/HoltLab/data/ica/012_ica.fif ...
Now restoring ICA solution ...
Ready.
# Running auto detection...
# Finding ICs that match patterns in EOG channels...
Using channel EXG3 as EOG channel
... filtering ICA sources
Setting up band-pass filter from 1 - 10 Hz
FIR filter parameters
---------------------
Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter:
- Windowed frequency-domain design (firwin2) method
- Hann window
- Lower passband edge: 1.00
- Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz)
- Upper passband edge: 10.00 Hz
- Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz)
- Filter length: 2048 samples (16.000 sec)
... filtering target
Setting up band-pass filter from 1 - 10 Hz
FIR filter parameters
---------------------
Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter:
- Windowed frequency-domain design (firwin2) method
- Hann window
- Lower passband edge: 1.00
- Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz)
- Upper passband edge: 10.00 Hz
- Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz)
- Filter length: 2048 samples (16.000 sec)
Using channel EXG4 as EOG channel ... filtering ICA sources Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec) ... filtering target Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec)
# Exclude EOG: {1, 34}
# Finding ICs that match patterns in ECG channels...
# Exclude ECG: set()
# Auto Exclude: {1, 34}
# Exclude list: [1, 34, 3, 4, 5, 7, 10, 25]
# Visualizing noise reduction...
Transforming to ICA space (35 components)
Zeroing out 8 ICA components
# Saving the updated data... Transforming to ICA space (35 components) Zeroing out 8 ICA components Overwriting existing file. Writing F:\HoltLab\data\raw_after_ica\012_after_ica.fif Closing F:\HoltLab\data\raw_after_ica\012_after_ica.fif [done] Done!
subject = '015'
ECG_list = []
noise_list = [5, 13, 17, 18, 29] # eye: 17, 29; muscle: 18; noise: 5, 13
reject(subject, EOG_channels, ECG_list, noise_list)
##### Processing subject 015 ...
# Reading Raw and ICA data...
Opening raw data file F:/HoltLab/data/filtered_raw/015_filtered_raw.fif...
Isotrak not found
Range : 0 ... 632447 = 0.000 ... 4940.992 secs
Ready.
Reading 0 ... 632447 = 0.000 ... 4940.992 secs...
Reading F:/HoltLab/data/ica/015_ica.fif ...
Now restoring ICA solution ...
Ready.
# Running auto detection...
# Finding ICs that match patterns in EOG channels...
Using channel EXG3 as EOG channel
... filtering ICA sources
Setting up band-pass filter from 1 - 10 Hz
FIR filter parameters
---------------------
Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter:
- Windowed frequency-domain design (firwin2) method
- Hann window
- Lower passband edge: 1.00
- Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz)
- Upper passband edge: 10.00 Hz
- Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz)
- Filter length: 2048 samples (16.000 sec)
... filtering target
Setting up band-pass filter from 1 - 10 Hz
FIR filter parameters
---------------------
Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter:
- Windowed frequency-domain design (firwin2) method
- Hann window
- Lower passband edge: 1.00
- Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz)
- Upper passband edge: 10.00 Hz
- Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz)
- Filter length: 2048 samples (16.000 sec)
Using channel EXG4 as EOG channel ... filtering ICA sources Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec) ... filtering target Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec)
Using channel EXG5 as EOG channel ... filtering ICA sources Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec) ... filtering target Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec)
# Exclude EOG: {0, 4, 29}
# Finding ICs that match patterns in ECG channels...
# Exclude ECG: set()
# Auto Exclude: {0, 4, 29}
# Exclude list: [0, 4, 5, 13, 17, 18, 29]
# Visualizing noise reduction...
Transforming to ICA space (37 components)
Zeroing out 7 ICA components
# Saving the updated data... Transforming to ICA space (37 components) Zeroing out 7 ICA components Overwriting existing file. Writing F:\HoltLab\data\raw_after_ica\015_after_ica.fif Closing F:\HoltLab\data\raw_after_ica\015_after_ica.fif [done] Done!
subject = '016'
ECG_list = []
noise_list = [1, 2, 4, 5, 7, 8] # eye: 7, 8; noise: 1, 2, 3, 5;
reject(subject, EOG_channels, ECG_list, noise_list)
##### Processing subject 016 ...
# Reading Raw and ICA data...
Opening raw data file F:/HoltLab/data/filtered_raw/016_filtered_raw.fif...
Isotrak not found
Range : 0 ... 658815 = 0.000 ... 5146.992 secs
Ready.
Reading 0 ... 658815 = 0.000 ... 5146.992 secs...
Reading F:/HoltLab/data/ica/016_ica.fif ...
Now restoring ICA solution ...
Ready.
# Running auto detection...
# Finding ICs that match patterns in EOG channels...
Using channel EXG3 as EOG channel
... filtering ICA sources
Setting up band-pass filter from 1 - 10 Hz
FIR filter parameters
---------------------
Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter:
- Windowed frequency-domain design (firwin2) method
- Hann window
- Lower passband edge: 1.00
- Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz)
- Upper passband edge: 10.00 Hz
- Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz)
- Filter length: 2048 samples (16.000 sec)
... filtering target
Setting up band-pass filter from 1 - 10 Hz
FIR filter parameters
---------------------
Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter:
- Windowed frequency-domain design (firwin2) method
- Hann window
- Lower passband edge: 1.00
- Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz)
- Upper passband edge: 10.00 Hz
- Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz)
- Filter length: 2048 samples (16.000 sec)
Using channel EXG4 as EOG channel ... filtering ICA sources Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec) ... filtering target Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec)
Using channel EXG5 as EOG channel ... filtering ICA sources Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec) ... filtering target Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec)
# Exclude EOG: {0, 4, 6}
# Finding ICs that match patterns in ECG channels...
# Exclude ECG: set()
# Auto Exclude: {0, 4, 6}
# Exclude list: [0, 1, 2, 4, 5, 6, 7, 8]
# Visualizing noise reduction...
Transforming to ICA space (35 components)
Zeroing out 8 ICA components
# Saving the updated data... Transforming to ICA space (35 components) Zeroing out 8 ICA components Overwriting existing file. Writing F:\HoltLab\data\raw_after_ica\016_after_ica.fif Closing F:\HoltLab\data\raw_after_ica\016_after_ica.fif [done] Done!
subject = '021'
ECG_list = []
noise_list = [2, 6, 11, 13, 30] # eye: 2; muscle: 6, 11, 13, 30
reject(subject, EOG_channels, ECG_list, noise_list)
##### Processing subject 021 ...
# Reading Raw and ICA data...
Opening raw data file F:/HoltLab/data/filtered_raw/021_filtered_raw.fif...
Isotrak not found
Range : 0 ... 756863 = 0.000 ... 5912.992 secs
Ready.
Reading 0 ... 756863 = 0.000 ... 5912.992 secs...
Reading F:/HoltLab/data/ica/021_ica.fif ...
Now restoring ICA solution ...
Ready.
# Running auto detection...
# Finding ICs that match patterns in EOG channels...
Using channel EXG3 as EOG channel
... filtering ICA sources
Setting up band-pass filter from 1 - 10 Hz
FIR filter parameters
---------------------
Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter:
- Windowed frequency-domain design (firwin2) method
- Hann window
- Lower passband edge: 1.00
- Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz)
- Upper passband edge: 10.00 Hz
- Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz)
- Filter length: 2048 samples (16.000 sec)
... filtering target
Setting up band-pass filter from 1 - 10 Hz
FIR filter parameters
---------------------
Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter:
- Windowed frequency-domain design (firwin2) method
- Hann window
- Lower passband edge: 1.00
- Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz)
- Upper passband edge: 10.00 Hz
- Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz)
- Filter length: 2048 samples (16.000 sec)
Using channel EXG4 as EOG channel ... filtering ICA sources Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec) ... filtering target Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec)
Using channel EXG5 as EOG channel ... filtering ICA sources Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec) ... filtering target Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec)
# Exclude EOG: set() # Finding ICs that match patterns in ECG channels... # Exclude ECG: set() # Auto Exclude: set() # Exclude list: [2, 6, 11, 13, 30] # Visualizing noise reduction... Transforming to ICA space (36 components) Zeroing out 5 ICA components
# Saving the updated data... Transforming to ICA space (36 components) Zeroing out 5 ICA components Overwriting existing file. Writing F:\HoltLab\data\raw_after_ica\021_after_ica.fif Closing F:\HoltLab\data\raw_after_ica\021_after_ica.fif [done] Done!
subject = '022'
ECG_list = []
noise_list = [0, 1, 2, 4, 7, 13, 14, 26, 29] # muscle: 14, 26, 29; noise: 0, 1, 2, 4, 7, 13
reject(subject, EOG_channels, ECG_list, noise_list)
##### Processing subject 022 ...
# Reading Raw and ICA data...
Opening raw data file F:/HoltLab/data/filtered_raw/022_filtered_raw.fif...
Isotrak not found
Range : 0 ... 660863 = 0.000 ... 5162.992 secs
Ready.
Reading 0 ... 660863 = 0.000 ... 5162.992 secs...
Reading F:/HoltLab/data/ica/022_ica.fif ...
Now restoring ICA solution ...
Ready.
# Running auto detection...
# Finding ICs that match patterns in EOG channels...
Using channel EXG3 as EOG channel
... filtering ICA sources
Setting up band-pass filter from 1 - 10 Hz
FIR filter parameters
---------------------
Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter:
- Windowed frequency-domain design (firwin2) method
- Hann window
- Lower passband edge: 1.00
- Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz)
- Upper passband edge: 10.00 Hz
- Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz)
- Filter length: 2048 samples (16.000 sec)
... filtering target
Setting up band-pass filter from 1 - 10 Hz
FIR filter parameters
---------------------
Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter:
- Windowed frequency-domain design (firwin2) method
- Hann window
- Lower passband edge: 1.00
- Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz)
- Upper passband edge: 10.00 Hz
- Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz)
- Filter length: 2048 samples (16.000 sec)
Using channel EXG4 as EOG channel ... filtering ICA sources Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec) ... filtering target Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec)
Using channel EXG5 as EOG channel ... filtering ICA sources Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec) ... filtering target Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec)
# Exclude EOG: {20}
# Finding ICs that match patterns in ECG channels...
# Exclude ECG: set()
# Auto Exclude: {20}
# Exclude list: [0, 1, 2, 4, 7, 13, 14, 20, 26, 29]
# Visualizing noise reduction...
Transforming to ICA space (37 components)
Zeroing out 10 ICA components
# Saving the updated data... Transforming to ICA space (37 components) Zeroing out 10 ICA components Overwriting existing file. Writing F:\HoltLab\data\raw_after_ica\022_after_ica.fif Closing F:\HoltLab\data\raw_after_ica\022_after_ica.fif [done] Done!
subject = '024'
ECG_list = []
noise_list = [30] # muscle: 30
reject(subject, EOG_channels, ECG_list, noise_list)
##### Processing subject 024 ...
# Reading Raw and ICA data...
Opening raw data file F:/HoltLab/data/filtered_raw/024_filtered_raw.fif...
Isotrak not found
Range : 0 ... 629503 = 0.000 ... 4917.992 secs
Ready.
Reading 0 ... 629503 = 0.000 ... 4917.992 secs...
Reading F:/HoltLab/data/ica/024_ica.fif ...
Now restoring ICA solution ...
Ready.
# Running auto detection...
# Finding ICs that match patterns in EOG channels...
Using channel EXG3 as EOG channel
... filtering ICA sources
Setting up band-pass filter from 1 - 10 Hz
FIR filter parameters
---------------------
Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter:
- Windowed frequency-domain design (firwin2) method
- Hann window
- Lower passband edge: 1.00
- Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz)
- Upper passband edge: 10.00 Hz
- Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz)
- Filter length: 2048 samples (16.000 sec)
... filtering target
Setting up band-pass filter from 1 - 10 Hz
FIR filter parameters
---------------------
Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter:
- Windowed frequency-domain design (firwin2) method
- Hann window
- Lower passband edge: 1.00
- Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz)
- Upper passband edge: 10.00 Hz
- Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz)
- Filter length: 2048 samples (16.000 sec)
Using channel EXG4 as EOG channel ... filtering ICA sources Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec) ... filtering target Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec)
Using channel EXG5 as EOG channel ... filtering ICA sources Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec) ... filtering target Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec)
# Exclude EOG: {1, 2}
# Finding ICs that match patterns in ECG channels...
# Exclude ECG: set()
# Auto Exclude: {1, 2}
# Exclude list: [1, 2, 30]
# Visualizing noise reduction...
Transforming to ICA space (35 components)
Zeroing out 3 ICA components
# Saving the updated data... Transforming to ICA space (35 components) Zeroing out 3 ICA components Overwriting existing file. Writing F:\HoltLab\data\raw_after_ica\024_after_ica.fif Closing F:\HoltLab\data\raw_after_ica\024_after_ica.fif [done] Done!
subject = '025'
ECG_list = []
noise_list = [0, 1, 4, 12, 28, 34] # eye: 12; muscle: 28; noise: 0, 1, 4, 34
reject(subject, EOG_channels, ECG_list, noise_list)
##### Processing subject 025 ...
# Reading Raw and ICA data...
Opening raw data file F:/HoltLab/data/filtered_raw/025_filtered_raw.fif...
Isotrak not found
Range : 0 ... 647679 = 0.000 ... 5059.992 secs
Ready.
Reading 0 ... 647679 = 0.000 ... 5059.992 secs...
Reading F:/HoltLab/data/ica/025_ica.fif ...
Now restoring ICA solution ...
Ready.
# Running auto detection...
# Finding ICs that match patterns in EOG channels...
Using channel EXG3 as EOG channel
... filtering ICA sources
Setting up band-pass filter from 1 - 10 Hz
FIR filter parameters
---------------------
Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter:
- Windowed frequency-domain design (firwin2) method
- Hann window
- Lower passband edge: 1.00
- Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz)
- Upper passband edge: 10.00 Hz
- Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz)
- Filter length: 2048 samples (16.000 sec)
... filtering target
Setting up band-pass filter from 1 - 10 Hz
FIR filter parameters
---------------------
Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter:
- Windowed frequency-domain design (firwin2) method
- Hann window
- Lower passband edge: 1.00
- Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz)
- Upper passband edge: 10.00 Hz
- Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz)
- Filter length: 2048 samples (16.000 sec)
Using channel EXG4 as EOG channel ... filtering ICA sources Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec) ... filtering target Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec)
Using channel EXG5 as EOG channel ... filtering ICA sources Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec) ... filtering target Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec)
# Exclude EOG: {2, 34}
# Finding ICs that match patterns in ECG channels...
# Exclude ECG: set()
# Auto Exclude: {2, 34}
# Exclude list: [0, 1, 2, 34, 4, 12, 28]
# Visualizing noise reduction...
Transforming to ICA space (37 components)
Zeroing out 7 ICA components
# Saving the updated data... Transforming to ICA space (37 components) Zeroing out 7 ICA components Overwriting existing file. Writing F:\HoltLab\data\raw_after_ica\025_after_ica.fif Closing F:\HoltLab\data\raw_after_ica\025_after_ica.fif [done] Done!
subject = '027'
ECG_list = []
noise_list = [0, 7, 8, 11, 12, 14, 26, 30, 36] # eye: 0; muscle: 7, 8, 11, 12, 14, 26, 30; noise: 36
reject(subject, EOG_channels, ECG_list, noise_list)
##### Processing subject 027 ...
# Reading Raw and ICA data...
Opening raw data file F:/HoltLab/data/filtered_raw/027_filtered_raw.fif...
Isotrak not found
Range : 0 ... 607871 = 0.000 ... 4748.992 secs
Ready.
Reading 0 ... 607871 = 0.000 ... 4748.992 secs...
Reading F:/HoltLab/data/ica/027_ica.fif ...
Now restoring ICA solution ...
Ready.
# Running auto detection...
# Finding ICs that match patterns in EOG channels...
Using channel EXG3 as EOG channel
... filtering ICA sources
Setting up band-pass filter from 1 - 10 Hz
FIR filter parameters
---------------------
Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter:
- Windowed frequency-domain design (firwin2) method
- Hann window
- Lower passband edge: 1.00
- Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz)
- Upper passband edge: 10.00 Hz
- Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz)
- Filter length: 2048 samples (16.000 sec)
... filtering target
Setting up band-pass filter from 1 - 10 Hz
FIR filter parameters
---------------------
Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter:
- Windowed frequency-domain design (firwin2) method
- Hann window
- Lower passband edge: 1.00
- Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz)
- Upper passband edge: 10.00 Hz
- Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz)
- Filter length: 2048 samples (16.000 sec)
Using channel EXG4 as EOG channel ... filtering ICA sources Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec) ... filtering target Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec)
Using channel EXG5 as EOG channel ... filtering ICA sources Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec) ... filtering target Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec)
# Exclude EOG: {8}
# Finding ICs that match patterns in ECG channels...
# Exclude ECG: set()
# Auto Exclude: {8}
# Exclude list: [0, 36, 7, 8, 11, 12, 14, 26, 30]
# Visualizing noise reduction...
Transforming to ICA space (37 components)
Zeroing out 9 ICA components
# Saving the updated data... Transforming to ICA space (37 components) Zeroing out 9 ICA components Overwriting existing file. Writing F:\HoltLab\data\raw_after_ica\027_after_ica.fif Closing F:\HoltLab\data\raw_after_ica\027_after_ica.fif [done] Done!
subject = '029'
ECG_list = []
noise_list = [0, 1, 3, 8, 13, 31, 32, 35] # eye: 0, 1, 8; muscle: 13, 32; noise: 3, 31, 35
reject(subject, EOG_channels, ECG_list, noise_list)
##### Processing subject 029 ...
# Reading Raw and ICA data...
Opening raw data file F:/HoltLab/data/filtered_raw/029_filtered_raw.fif...
Isotrak not found
Range : 0 ... 675455 = 0.000 ... 5276.992 secs
Ready.
Reading 0 ... 675455 = 0.000 ... 5276.992 secs...
Reading F:/HoltLab/data/ica/029_ica.fif ...
Now restoring ICA solution ...
Ready.
# Running auto detection...
# Finding ICs that match patterns in EOG channels...
Using channel EXG3 as EOG channel
... filtering ICA sources
Setting up band-pass filter from 1 - 10 Hz
FIR filter parameters
---------------------
Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter:
- Windowed frequency-domain design (firwin2) method
- Hann window
- Lower passband edge: 1.00
- Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz)
- Upper passband edge: 10.00 Hz
- Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz)
- Filter length: 2048 samples (16.000 sec)
... filtering target
Setting up band-pass filter from 1 - 10 Hz
FIR filter parameters
---------------------
Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter:
- Windowed frequency-domain design (firwin2) method
- Hann window
- Lower passband edge: 1.00
- Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz)
- Upper passband edge: 10.00 Hz
- Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz)
- Filter length: 2048 samples (16.000 sec)
Using channel EXG4 as EOG channel ... filtering ICA sources Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec) ... filtering target Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec)
Using channel EXG5 as EOG channel ... filtering ICA sources Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec) ... filtering target Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec)
# Exclude EOG: {8, 7}
# Finding ICs that match patterns in ECG channels...
# Exclude ECG: set()
# Auto Exclude: {8, 7}
# Exclude list: [0, 1, 32, 3, 35, 7, 8, 13, 31]
# Visualizing noise reduction...
Transforming to ICA space (36 components)
Zeroing out 9 ICA components
# Saving the updated data... Transforming to ICA space (36 components) Zeroing out 9 ICA components Overwriting existing file. Writing F:\HoltLab\data\raw_after_ica\029_after_ica.fif Closing F:\HoltLab\data\raw_after_ica\029_after_ica.fif [done] Done!
subject = '030'
ECG_list = []
noise_list = [0, 3, 9, 13, 22, 28, 35] # eye: 0, 3, 22; muscle: 9, 13; noise: 28, 35
reject(subject, EOG_channels, ECG_list, noise_list)
##### Processing subject 030 ...
# Reading Raw and ICA data...
Opening raw data file F:/HoltLab/data/filtered_raw/030_filtered_raw.fif...
Isotrak not found
Range : 0 ... 671359 = 0.000 ... 5244.992 secs
Ready.
Reading 0 ... 671359 = 0.000 ... 5244.992 secs...
Reading F:/HoltLab/data/ica/030_ica.fif ...
Now restoring ICA solution ...
Ready.
# Running auto detection...
# Finding ICs that match patterns in EOG channels...
Using channel EXG3 as EOG channel
... filtering ICA sources
Setting up band-pass filter from 1 - 10 Hz
FIR filter parameters
---------------------
Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter:
- Windowed frequency-domain design (firwin2) method
- Hann window
- Lower passband edge: 1.00
- Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz)
- Upper passband edge: 10.00 Hz
- Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz)
- Filter length: 2048 samples (16.000 sec)
... filtering target
Setting up band-pass filter from 1 - 10 Hz
FIR filter parameters
---------------------
Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter:
- Windowed frequency-domain design (firwin2) method
- Hann window
- Lower passband edge: 1.00
- Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz)
- Upper passband edge: 10.00 Hz
- Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz)
- Filter length: 2048 samples (16.000 sec)
Using channel EXG4 as EOG channel ... filtering ICA sources Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec) ... filtering target Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec)
Using channel EXG5 as EOG channel ... filtering ICA sources Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec) ... filtering target Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec)
# Exclude EOG: {1, 2}
# Finding ICs that match patterns in ECG channels...
# Exclude ECG: set()
# Auto Exclude: {1, 2}
# Exclude list: [0, 1, 2, 3, 35, 9, 13, 22, 28]
# Visualizing noise reduction...
Transforming to ICA space (36 components)
Zeroing out 9 ICA components
# Saving the updated data... Transforming to ICA space (36 components) Zeroing out 9 ICA components Overwriting existing file. Writing F:\HoltLab\data\raw_after_ica\030_after_ica.fif Closing F:\HoltLab\data\raw_after_ica\030_after_ica.fif [done] Done!
subject = '031'
ECG_list = []
noise_list = [0, 4, 5, 6, 7, 12, 16, 35] # eye: 0, 5; muscle: 7, 12, 16; noise: 4, 6, 35
reject(subject, EOG_channels, ECG_list, noise_list)
##### Processing subject 031 ...
# Reading Raw and ICA data...
Opening raw data file F:/HoltLab/data/filtered_raw/031_filtered_raw.fif...
Isotrak not found
Range : 0 ... 777471 = 0.000 ... 6073.992 secs
Ready.
Reading 0 ... 777471 = 0.000 ... 6073.992 secs...
Reading F:/HoltLab/data/ica/031_ica.fif ...
Now restoring ICA solution ...
Ready.
# Running auto detection...
# Finding ICs that match patterns in EOG channels...
Using channel EXG3 as EOG channel
... filtering ICA sources
Setting up band-pass filter from 1 - 10 Hz
FIR filter parameters
---------------------
Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter:
- Windowed frequency-domain design (firwin2) method
- Hann window
- Lower passband edge: 1.00
- Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz)
- Upper passband edge: 10.00 Hz
- Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz)
- Filter length: 2048 samples (16.000 sec)
... filtering target
Setting up band-pass filter from 1 - 10 Hz
FIR filter parameters
---------------------
Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter:
- Windowed frequency-domain design (firwin2) method
- Hann window
- Lower passband edge: 1.00
- Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz)
- Upper passband edge: 10.00 Hz
- Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz)
- Filter length: 2048 samples (16.000 sec)
Using channel EXG4 as EOG channel ... filtering ICA sources Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec) ... filtering target Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec)
Using channel EXG5 as EOG channel ... filtering ICA sources Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec) ... filtering target Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec)
# Exclude EOG: {1, 2}
# Finding ICs that match patterns in ECG channels...
# Exclude ECG: set()
# Auto Exclude: {1, 2}
# Exclude list: [0, 1, 2, 35, 4, 5, 6, 7, 12, 16]
# Visualizing noise reduction...
Transforming to ICA space (36 components)
Zeroing out 10 ICA components
# Saving the updated data... Transforming to ICA space (36 components) Zeroing out 10 ICA components Overwriting existing file. Writing F:\HoltLab\data\raw_after_ica\031_after_ica.fif Closing F:\HoltLab\data\raw_after_ica\031_after_ica.fif [done] Done!
subject = '032'
ECG_list = []
noise_list = [0, 1, 2, 18, 19, 20, 24, 30] # eye: 24; muscle: 18, 19, 20; noise: 0, 1, 2, 14, 30;
reject(subject, EOG_channels, ECG_list, noise_list)
##### Processing subject 032 ...
# Reading Raw and ICA data...
Opening raw data file F:/HoltLab/data/filtered_raw/032_filtered_raw.fif...
Isotrak not found
Range : 0 ... 663551 = 0.000 ... 5183.992 secs
Ready.
Reading 0 ... 663551 = 0.000 ... 5183.992 secs...
Reading F:/HoltLab/data/ica/032_ica.fif ...
Now restoring ICA solution ...
Ready.
# Running auto detection...
# Finding ICs that match patterns in EOG channels...
Using channel EXG3 as EOG channel
... filtering ICA sources
Setting up band-pass filter from 1 - 10 Hz
FIR filter parameters
---------------------
Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter:
- Windowed frequency-domain design (firwin2) method
- Hann window
- Lower passband edge: 1.00
- Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz)
- Upper passband edge: 10.00 Hz
- Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz)
- Filter length: 2048 samples (16.000 sec)
... filtering target
Setting up band-pass filter from 1 - 10 Hz
FIR filter parameters
---------------------
Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter:
- Windowed frequency-domain design (firwin2) method
- Hann window
- Lower passband edge: 1.00
- Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz)
- Upper passband edge: 10.00 Hz
- Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz)
- Filter length: 2048 samples (16.000 sec)
Using channel EXG4 as EOG channel ... filtering ICA sources Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec) ... filtering target Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec)
Using channel EXG5 as EOG channel ... filtering ICA sources Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec) ... filtering target Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec)
# Exclude EOG: {8, 16, 6}
# Finding ICs that match patterns in ECG channels...
# Exclude ECG: set()
# Auto Exclude: {8, 16, 6}
# Exclude list: [0, 1, 2, 6, 8, 16, 18, 19, 20, 24, 30]
# Visualizing noise reduction...
Transforming to ICA space (34 components)
Zeroing out 11 ICA components
# Saving the updated data... Transforming to ICA space (34 components) Zeroing out 11 ICA components Overwriting existing file. Writing F:\HoltLab\data\raw_after_ica\032_after_ica.fif Closing F:\HoltLab\data\raw_after_ica\032_after_ica.fif [done] Done!
subject = '033'
ECG_list = []
noise_list = [1, 2, 3, 4, 6, 7, 12, 14, 18, 20, 28, 31, 35] # eye: 2, 3, 6, 7, 12, 20, 31, 35; muscle: 1, 4, 18, 28; noise: 14
reject(subject, EOG_channels, ECG_list, noise_list)
##### Processing subject 033 ...
# Reading Raw and ICA data...
Opening raw data file F:/HoltLab/data/filtered_raw/033_filtered_raw.fif...
Isotrak not found
Range : 0 ... 650495 = 0.000 ... 5081.992 secs
Ready.
Reading 0 ... 650495 = 0.000 ... 5081.992 secs...
Reading F:/HoltLab/data/ica/033_ica.fif ...
Now restoring ICA solution ...
Ready.
# Running auto detection...
# Finding ICs that match patterns in EOG channels...
Using channel EXG3 as EOG channel
... filtering ICA sources
Setting up band-pass filter from 1 - 10 Hz
FIR filter parameters
---------------------
Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter:
- Windowed frequency-domain design (firwin2) method
- Hann window
- Lower passband edge: 1.00
- Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz)
- Upper passband edge: 10.00 Hz
- Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz)
- Filter length: 2048 samples (16.000 sec)
... filtering target
Setting up band-pass filter from 1 - 10 Hz
FIR filter parameters
---------------------
Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter:
- Windowed frequency-domain design (firwin2) method
- Hann window
- Lower passband edge: 1.00
- Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz)
- Upper passband edge: 10.00 Hz
- Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz)
- Filter length: 2048 samples (16.000 sec)
Using channel EXG4 as EOG channel ... filtering ICA sources Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec) ... filtering target Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec)
Using channel EXG5 as EOG channel ... filtering ICA sources Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec) ... filtering target Setting up band-pass filter from 1 - 10 Hz FIR filter parameters --------------------- Designing a two-pass forward and reverse, zero-phase, non-causal bandpass filter: - Windowed frequency-domain design (firwin2) method - Hann window - Lower passband edge: 1.00 - Lower transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 0.75 Hz) - Upper passband edge: 10.00 Hz - Upper transition bandwidth: 0.50 Hz (-12 dB cutoff frequency: 10.25 Hz) - Filter length: 2048 samples (16.000 sec)
# Exclude EOG: {9, 6}
# Finding ICs that match patterns in ECG channels...
# Exclude ECG: set()
# Auto Exclude: {9, 6}
# Exclude list: [1, 2, 3, 4, 35, 6, 7, 9, 12, 14, 18, 20, 28, 31]
# Visualizing noise reduction...
Transforming to ICA space (36 components)
Zeroing out 14 ICA components
# Saving the updated data... Transforming to ICA space (36 components) Zeroing out 14 ICA components Overwriting existing file. Writing F:\HoltLab\data\raw_after_ica\033_after_ica.fif Closing F:\HoltLab\data\raw_after_ica\033_after_ica.fif [done] Done!